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Abstract 

Results are reported from a search for the production of three-jet resonances in pp 
collisions at a center-of-mass energy \/s = 7TeV. The study uses the data sample 
collected by the CMS experiment at the LHC in 2011, corresponding to an integrated 
luminosity of 5.0 fb -1 . Events with high jet multiplicity and a large scalar sum of 
jet transverse momenta are analyzed for the presence of resonances in the three-jet 
invariant mass spectrum. The number of events observed is found to be in agree- 
ment with the expectation from standard model processes. Limits are set on the cross 
section for gluino pair prediction in an R-parity-violating supersymmetry model, for 
gluino masses greater than 280 GeV. Assuming a branching fraction for gluino decay 
into three jets of 100%, gluino masses below 460 GeV are excluded at 95% confidence 
level. These results significantly extend the range of previous limits. 
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Many extensions of the standard model (SM) predict the existence of new, strongly interact- 
ing particles that lead to final states with high jet multiplicities. Examples of such particles 
include gluinos in R-parity-violating (RPV) SUSY models d [3 and heavy colored fermion 
resonances, which transform as octets under SU(3)c USHSI- Despite large production cross sec- 
tions predicted by many of these models, the vast majority of current experimental searches 
are insensitive to them because they use signal selection criteria based on low jet multiplicities, 
such as dijet resonances, or require the presence of large missing transverse energy. While this 
strategy has proved effective in establishing cross section limits for many models of physics be- 
yond the SM up to TeV mass scales EL the currently-published limits for models predicting 
high jet multiplicity final states |9l-[TTf are much less stringent. Specifically, there are only two 
published results from searches involving six jet final states at hadron colliders, one from the 
Tevatron [12] and the other from the Large Hadron Collider (LHC) [13J, which exclude RPV 
gluinos with masses below 144 GeV and between 200 and 280 GeV, respectively. 

This Letter presents the results of a search for three-jet resonances in multijet events in proton- 
proton (pp) collisions at a/s = 7 TeV. The event sample was collected using the Compact Muon 
Solenoid (CMS) detector [14J at the LHC. The integrated luminosity is 4.98 ± 0.11 fb _1 llTBI . 
corresponding to the CMS data sample recorded in 2011. In this extension of previous CMS 
results ||T3| using 35 pb _1 of data, events with at least six high-transverse-momentum (px) jets 
are investigated for evidence of three-jet resonances using the jet-ensemble technique ||T2lH3l . 

The CMS detector is a multi-purpose apparatus described in detail in Ref. |(T4| . A high-resolution 
silicon pixel and strip tracker, immersed in the 3.8 T magnetic field of the superconducting 
solenoid, provides charged particle tracking coverage for \tj\ < 2.4, where rj = — ln[tan(0/2)] 
is the pseudorapidity and 9 is the polar angle measured with respect to the counterclockwise 
proton beam direction. Jet energy deposits are measured using electromagnetic (ECAL) and 
hadronic (HCAL) calorimeters. The finely-segmented lead-tungstate crystal ECAL consists of 
a barrel and two endcap regions. The ECAL barrel covers the pseudorapidity range \tj\ < 1.4 
with a granularity of At/ x A<p = 0.0174 x 0.0174, where cp is the azimuthal angle, measured in 
radians, and the endcaps cover 1.4 < \rj\ < 3.0. The HCAL extends out to an \rj \ of about 5.0. 
Its central and endcap regions consist of brass/ scintillator sampling calorimeters that cover 
\t] | < 3.0 with a granularity of At] x Acp = 0.087 x 0.087 for central rapidities. Further cov- 
erage is provided by a steel /quartz-fiber Cherenkov calorimeter. Muons are measured in gas 
detectors embedded in the steel return yoke outside the solenoid. 

Events are recorded using a two-tiered trigger system [14j. Objects satisfying the requirements 
at the first level are passed to the high-level trigger (HLT) which further reduces the total 
recorded rate. Triggers based on the scalar sum of all transverse energy from jets (Hj), re- 
constructed using only calorimeter information, are used to select events. To reduce the effects 
of multiple pp interactions in the same bunch crossings (pileup), which produce a large number 
of low-px particles, jets considered in the HLT selection are required to have pj > 40 GeV. The 
HLT threshold for Hj ranged from 350 to 750 GeV over the course of the 2011 data collection 
period, depending on the instantaneous luminosity. 

A further set of criteria is imposed on events passing the trigger. Events are required to contain 
at least one well-reconstructed primary event vertex [[161 . Jets are required to have pj > 70 GeV 
and 1 1] | < 3.0, which suppresses background and minimizes pileup effects. Events with pair- 
produced gluinos are expected to have high jet multiplicity and large values of Hj. We there- 
fore require events to contain at least six jets, and to have a total scalar sum of jet pj exceeding 
900 GeV. The latter requirement ensures that the trigger is fully efficient for these events. 

Individual objects reconstructed using the CMS particle-flow algorithm [17J serve as input for 
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jet reconstruction. The particle-flow algorithm |17) combines calorimeter information with re- 
constructed tracks to identify individual particles such as photons, leptons, and both neutral 
and charged hadrons within jets. The energy of photons is obtained directly from the calibrated 
ECAL measurement. The determination of the energy of electrons comes from a combination 
of the track momentum at the main interaction vertex, the corresponding ECAL cluster energy, 
and the sum of all energy from bremsstrahlung photons associated to the track. The energy of 
muons is obtained from the corresponding track momentum. Charged hadron energy is deter- 
mined from a combination of the track momentum and the corresponding ECAL and HCAL 
energies, corrected for hardware zero-suppression effects, and calibrated for the non-linear re- 
sponse of the HCAL. Finally, the energy of neutral hadrons is obtained from the corresponding 
calibrated ECAL and HCAL energies. 

Jet reconstruction is performed using the anti-fcx algorithm [18] with a distance parameter of 
0.5 in rj-(p space. Jet energy scale corrections [19] derived from data and Monte Carlo (MC) 
simulation are applied to account for the non-linear and non-uniform response of the HCAL. 
A small residual correction factor is included to correct for differences in jet response between 
the detector and its simulation. The combined corrections are on the order of 5-10%, and their 
corresponding uncertainties range from 1 to 6%, depending on the jet pseudorapidity and px- 
Jet quality criteria 11201 are applied to remove misidentified jets, which arise primarily from 
calorimeter noise. For both data and simulated signal events, more than 99.8% of all selected 
jets satisfy these criteria. After imposing all selection requirements, 114599 events remain. 

The benchmark signal model for the analysis features pair-produced gluinos, whose produc- 
tion and decay are simulated using the pythia l|2"T1| v6.424 event simulation program with 
the CTEQ6L1 parton distribution function (PDF) l[22| . The gluino decay is modeled as an ef- 
fective RPV coupling between the gluino and three light quarks, with a branching fraction 
of 100%. The gluino mass is varied between 250 and 500 GeV in steps of 50 GeV, and addi- 
tional mass points are generated at 750, 1000, 1250, and 1500 GeV. The leading-order (LO) 
cross section from PYTHIA is 92 pb for a gluino mass of 250 GeV and 2.1 pb for a gluino mass 
of 450 GeV, falling to about 1.3 x 10~ 5 pb for a gluino mass of 1500 GeV. All superpartners ex- 
cept the gluino are decoupled, the natural width of the gluino resonance is taken to be much 
smaller than the resolution of the detector, and no intermediate particles are produced in the 
gluino decay. Next-to-leading-order (NLO) correction factors (K-factors) are calculated using 
the PROSPINO |23| program and are applied to the LO cross sections. These values rise from 1.8 
for a gluino mass of 250 GeV to 2.1 for a gluino mass of 450 GeV, and to 4.6 for a gluino mass of 
1500 GeV. Simulation of the CMS detector response is performed using Geant4 [24]. 

The principal background arises from multijet events produced in strong interaction (QCD) 
processes. All other SM processes, such as fully-hadronic decays of top quark pairs, are pre- 
dicted to make a negligible contribution. The QCD background arises primarily from events 
with two high-px jets from a hard interaction, combined with gluon jets from initial- or final- 
state radiation. We use the jet-ensemble technique to reconstruct the gluino candidates and 
to suppress background. The six highest-px jets in each event are combined in all possible 
three-jet combinations (triplets), resulting in 20 unique triplets per event. For signal events, 
at most two of the triplets correspond to a parent gluino, while the other triplets constitute a 
combinatorial background. Thus, background arises not only from QCD multijet events, but 
also from the signal events themselves. To improve sensitivity to the presence of a three-jet 
resonance, a requirement that exploits the observed linear correlation between the triplet mass 
and the scalar sum of the three jets' px values in background triplet combinations is imposed. 
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Figure 1: Simulated triplet jet mass Mjjj versus the triplet scalar pj of all 20 triplets for a gluino 
mass of 400 GeV, shown as a contoured distribution, with that of the triplets whose jets all 
originate from the same gluino parent overlaid. The expectation for the QCD background is 
shown in the inset. All triplets falling to the right of the dashed line pass the requirement of 
Eq. Q for a value of A = 160 GeV. 

We require 

M )jj<I>T-A, (1) 

i=l 

where M™ is the triplet mass, Y%=\ p' T is the sum of the magnitudes of the transverse momenta 
of the jets within the triplet, and A is an offset adjusted to optimize signal sensitivity. Note 
that for signal triplets, Mjjj is fixed by the gluino mass and thus independent of its three jets' 
scalar pj sum. Figure [T] shows the simulated triplet mass versus the triplet scalar pj for a 
gluino mass of 400 GeV. The value of A is determined by maximizing the ratio of the number 
of signal triplets to the sum of the number of signal and background triplets in a two-standard- 
deviation (a) interval about the mean of a Gaussian function that describes the gluino mass 
peak. A common value A = 160 GeV is used for all gluino masses considered, and typically at 
most one triplet per event survives the requirement from Eq. flj. The fraction of signal triplets, 
whose jets all originate from the same parent gluino, that pass this requirement ranges from 2 to 
13% depending on the gluino mass, while the efficiency to select all other triplet combinations 
in a signal event is below 1%. 

After the final selection, residual background remains from both QCD multijet events and the 
combinatorial background triplets in possible gluino signal events. The latter contribute only 
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minimally, and the shape of their mass distribution is found to be similar to that of the dom- 
inant QCD multijet background. We therefore consider the two background components to- 
gether and parameterize their shape using the smoothly falling distribution [25 j given by: 

da_ PoCl-AVv/s)^ 

The parameters P , Pi, P 2 , and P 3 have values of (61.0 ± 1.8) x 10 3 GeV _1 , 53.8 ± 0.7, -1.90 ± 
0.04, and —0.13 ± 0.01, respectively, as determined by fitting this function to the data sample in 
the range from 260 to 1625 GeV using a x 2 minimization technique. Thus, any potential signal 
would manifest itself as a localized positive deviation from the expectation of the background 
fit. 

To estimate the number of signal events expected after all selection criteria are applied, the sum 
of a Gaussian function that represents the signal and the four-parameter function (Eq. ([2])) that 
models the background are fitted to the simulated gluino M™ distribution, for Mj-j > 260 GeV. 
Below 260 GeV, the triplet mass distribution falls off rapidly due to restrictions imposed by 
the Hj threshold in the HLT The width of the Gaussian function modeling the signal varies 
according to the detector resolution from 17 GeV for a gluino mass of 250 GeV, to 100 GeV for 
a gluino mass of 1500 GeV, as determined from simulation. The integral of the Gaussian com- 
ponent provides the estimate for the expected number of signal triplets reconstructed, and the 
value of this integral, divided by the number of signal events generated, determines the signal 
acceptance for each gluino mass. The signal acceptance per event includes all selection criteria 
and is parameterized with a second-order polynomial as a function of gluino mass. The event 
acceptance is 0.25% for a gluino of mass 250 GeV, 1.5% for a gluino of mass 450 GeV, and rises 
to 2.6% for a gluino of mass 1500 GeV. 

The systematic uncertainty on the signal acceptance is evaluated as follows. The uncertainty 
related to the jet energy scale [19] is evaluated by varying the jet energy scale correction within 
its uncertainties, then recalculating the acceptance for different gluino mass values. The largest 
difference with respect to the nominal acceptance is taken as the systematic uncertainty and 
ranges from 9 to 18%. The difference between the calculated gluino acceptance at each mass 
point and the parameterized acceptance is taken as an uncertainty for each mass value and 
varies from 1 to 18%. The level of initial- and final-state radiation is varied following a prescrip- 
tion used in many other results 1126 1 , where the relative amount of each is separately increased 
and decreased with respect to the nominal value. The associated uncertainty is evaluated in a 
similar manner to that described for the jet energy scale uncertainties. The difference of 3 to 
5% with respect to the nominal acceptance is taken as the systematic uncertainty. To determine 
the effect of pileup on the signal acceptance, the MC signal samples are reweighted such that 
the distribution of reconstructed primary vertices is shifted high and then low by one standard 
deviation compared to the nominal analysis (in the nominal analysis, the MC distribution is 
reweighted to correspond to the measured distribution). The acceptance procedure is repeated 
at both points and the largest difference with respect to the nominal acceptance is taken as 
the systematic uncertainty, which varies from 1 to 3%. The uncertainty corresponding to the 
choice of PDF set is evaluated by reweighting the events to correspond to each of its associated 
eigenvector sets and repeating the acceptance calculation. The difference, added in quadrature 
from all sets, is 4% and is taken as the systematic uncertainty. These contributions, combined 
with those from the luminosity measurement (2.2%), yield a total systematic uncertainty on the 
signal acceptance between 14 and 26%, depending on the gluino mass. 

The measured mass distribution for the selected jet triplets is shown in Fig. |2j where the four- 
parameter background fit is represented by the solid line. The expectation for a 300 GeV and 
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Figure 2: Mass distribution for selected jet triplets. The four-parameter background fit is shown 
by the solid line. The NLO expectation for a 300 GeV gluino signal is shown by the dotted line, 
and that of a 450 GeV gluino signal by the dash-dotted line. Each of the two is normalized to 
the integrated luminosity of the data sample. 

a 450 GeV gluino signal are also included as a dotted and dash-dotted line, respectively, and 
each are normalized to the integrated luminosity of the data sample. We choose the width of 
bins to equal 10% of the mass value. This yields bin widths in general correspondence with the 
expected instrumental mass resolution. 

To search for an excess in the number of triplets with respect to the expected background, we 
compute the difference between the measured number of triplets observed in each mass bin 
and the number of triplets predicted in this bin by the background fit. Figure [3] shows this 
difference, normalized to the background prediction (Fig. 3a) and to the uncertainty in the 
background prediction (Fig. 3b). In both cases, the predicted background is in agreement with 
the data over the full mass range. 

Because no excess of events is observed in the measured mass distribution, we proceed to 
place upper limits on the cross section times the branching fraction for the production of three- 
jet resonances in the data. Both observed and expected limits are calculated using a modified 
frequentist CL S method with a one-sided profile likelihood test statistic [27J. The background 
model parameters and their corresponding uncertainties are taken from the fit of the four- 
parameter background function (Eq. ([2])) to the data. The uncertainties on the parameters that 
describe the background shape are included as log-normal constraints, and they are 1.3% for P\, 
2.2% for Pi, 8.6% for P3, with a normalization uncertainty on Pq of 3%. The central value is set 
to the best fit value and the width to one standard deviation. The range is truncated at ±3c. In 
addition to the background parameters, log-normal constraints are included for the acceptance 
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Figure 3: Difference between the measured triplet mass distribution and the fitted background 
parametrization, divided by the fitted value (a) or by the statistical uncertainty S on the fitted 
value (b), shown for both data and the NLO expectation of two gluino models, one with a 
gluino mass of 300 GeV and the other with a gluino mass of 450 GeV. 

and integrated luminosity. These systematic uncertainties are treated as global constraints in 
the likelihood. 

The observed and expected 95% confidence level (CL) upper limits on the gluino pair produc- 
tion cross section times branching fraction as a function of gluino mass are presented in Fig. |4] 
The corresponding 95% CL lower limit on the gluino mass is determined by finding the mass 
value at which the limit line crosses that of the NLO gluino cross section, assuming a branching 
fraction for the signal model of 100%. We perform the search in the region of M™ above 280 GeV 
and exclude masses below 460 GeV at 95% CL. 

In summary, a search has been performed for three-jet resonance production in pp collisions 
at a center-of-mass energy of 7 TeV, using a data sample corresponding to an integrated lumi- 
nosity of 5.0 fb . Events having the properties of high jet multiplicity and large scalar sum 
of jet pi, which are expected signatures of high-mass hadronic resonances, were analyzed for 
the presence of signal events. The number of observed events was found to be in agreement 
with the expectation from standard model processes. For a mass range above 280 GeV, the pro- 
duction of gluinos, modeled as an effective RPV coupling between the gluino and three light 
quarks with a branching fraction of 100%, has been excluded for masses below 460 GeV at 95% 
CL. These are the most stringent limits to date. 
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CMS 5.0 fb 1 = 7 TeV 
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Figure 4: Observed and expected 95% CL upper limits on the cross section times branching frac- 
tion for gluino pair production followed by RPV decay of each gluino to three light-flavored 
quark jets. Also shown are the ±lc and ±2cj bands on the expected limit, as well as the the- 
oretical LO and NLO cross sections for gluino production, assuming a branching fraction of a 
gluino decay into three jets of 100%. 
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